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LENTICULAR- PRINTING CALIBRATION TARGETS 

5 

Field of the Present Invention 

The present invention relates to lenticular printing, 
and more particular to lenticular printing calibration 
targets . 

10 

Background of the Present Invention 

Plastic-covered, ribbed, pictures that appear to move 
or change scenes are familiar to most Americans. These 
special effects can also include a three-dimensional 

15 effect. Lenticular sheets are laminated on to flat 

printed image. The cylindrical lenses are arranged in 
parallel strips and allow only one of the interlaced 
images to be seen by a viewer at any one angle. But as 
the angle of view changes, the scene can change and give 

20 the illusion of three-dimensions, morphing, animation, and 
flipping . 

The lenticular sheets typically include a series of 
optical grade cylindrical lenses, or lenticules. In use, 
several images can be printed underneath the lenticules. 

25 Each image is divided into strips and one strip from each 
image is printed underneath each lenticule. The 
lenticules refract the different strips so that one of the 
images is visible at a time according to the angle of 
observation. Changing the angle of observation causes 

30 different images to appear and in consequence the image 
appears to change . 
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For lenticular printing to be effective the image 
strips must be correctly aligned with the lenticules. 
Sheets of lenticular film are typically manufactured for a 
specified number of lines (or lens) per inch rating which 
5 can range 6-3 00 lines per inch dpi) . The pitch of the 
lenticules is the distance between lenticule centers. 
Generally, the pitch of the lenticule is a function of the 
desired viewing distance. A higher pitch is used for 
longer viewing distances. 

10 Effective lenticular printing requires extremely 

accurate pitch calibration of the lenticules, e.g., up to 
l/100 th of the lines-per-inch pitch. Due to imperfections 
in the manufacturing process, the actual pitch varies from 
the nominal value. Furthermore, it is important that the 

15 errors in the pitch of the interlaced print below the 
lenticules do not aggregate. 

Summary Of The Present Invention 

20 

Briefly, a method embodiment of the present invention 
generates lenticular-printing calibration targets. A 
square wave is generated at a frequency determined by a 
target pitch for the lenticules. The square wave is 
25 filtered to eliminate aliased harmonics related to a 
Nyquist frequency determined by the resolution of the 
printer and the required calibration precision. The 
resulting spatial domain square wave is printed for a 
lenticular sheet to produce a calibration target. 

30 
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Brief Description Of The Drawings 



Fig. 1 is a flowchart diagram of a method embodiment 
of the present invention for generating lenticular- 
5 printing calibration targets; 

Fig. 2A is a spatial domain plot of a 50% duty cycle 
square wave, the vertical axis is grayscale level between 
white and black; 

Fig. 2B is a frequency domain representation of the 
10 square wave of Fig. 2A; 

Fig. 3A is a frequency domain representation of Fig. 
2B after filtering; 

Fig. 3B is the spatial domain target plot of the 
waveform represented by Fig. 3A; 
15 Fig. 4A is a frequency domain representation of Fig. 

3B after eliminating undesirable oscillations by clamping; 

Fig. 4B is the spatial domain target plot of the 
final target waveform represented by Fig. 4A; and 

Fig. 5 is a cross section diagram of a lenticular 
20 printed image showing how a calibration target can be 

viewed from different angles by the right and left eyes. 

Detailed Description of the Embodiments 

25 

Fig. 1 represents a method embodiment of the present 
invention for generating lenticular-printing calibration 
targets, and is referred to herein by the general 
reference numeral 100. For example, the nominal pitch of 
30 a lenticular sheet is specified to be 20.0 lines-per-inch, 
the actual product will vary. Calibration targets can be 
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used to measure the true pitch, e.g., with calibration 
precision of l/10 th of an inch. 

Method 100 comprises generating a square wave in a 
step 101. A fast Fourier transform in a step 102 converts 
5 this to the frequency domain. A step 104 is then easily 
able to remove aliasing frequencies. A step 106 
reconverts to the time domain with an inverse fast Fourier 
transform. A step 108 clamps the resulting image. And a 
step 110 prints the calibration targets on sheets of 

10 paper. Users thereafter measure the actual pitch of a 

lenticular sheet by seeing which calibration target sheet 
produces the maximum effect. 

The software driving laser and ink jet printers can be 
used to scale the sizes of images output to a printer. 

15 Printer resolutions can well exceed 3 00 dots-per-inch 

(dpi). Calibration targets are printed that range 19.0- 
21.0 lenticular lines-per-inch, in increments of 0.1. The 
true pitch of a lenticular sheet can be determined by 
viewing the line patterns generated by each calibration 

20 target through the lenticular sheet, and observing which 
pattern switches best between all black and all white, 
over the range of viewing angles. 

A typical calibration target will be for 20.5 lines- 
per-inch, and a typical printer resolution will be 300 

25 dpi. A square wave generator with a 50% duty cycle, and 
an amplitude 1.0, is generated at a frequency dictated by 
a target pitch of 20.5 lines-per-inch. 

A step 102 involves computing the fast Fourier 
transform (FFT) of the generated square wave. The number 

30 of samples is determined by the resolution of the printer 
and the calibration precision required, that is the 
product 10 x 300 dpi, or 3000 samples. As illustrated in 
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Fig. 2B, the FFT 202 of this square wave produces a 
fundamental component peak 204 at sample 205. However it 
also produces odd harmonics 206 at odd multiples of the 
fundamental, that is at samples 615, 1025 and 1435. 
5 Aliasing 208 is also produced at neighboring frequencies. 

To avoid the generation of moire frequencies in the 
spatial domain, only the fundamental and the harmonics 
below the Nyquist frequency are required. The Nyquist 
frequency is the maximum frequency that is represented in 

10 the Fourier transform without aliasing,. Its location 
occurs at half the length of the FFT used. Thus the 
Nyquist frequency is 150 dpi represented by the 1500 th 
sample of the FFT. As a result, step 104 involves zeroing 
out such other frequencies. 

15 The result is a filtered transform which can be seen 

in Fig. 3A at 300. All other FFT components are due to 
aliased harmonics and are eliminated in this process. 

The conjugate frequency is computed and retained. 
The FFT samples at the corresponding conjugate locations 

20 of (3000-205), (3000-615), (3000-1025) and (3000-1435) are 
also retained to assure that the inverse transform is 
real . 

The inverse FFT is then computed in step 106 to 
reconstruct the square wave in the spatial domain, the 
25 imaginary component is discarded and the real part 
retained. As illustrated in Fig. 3B, this procedure 
results in an interpolated set of bars 302 that are close 
to being square waves but incorporate small oscillations 
304. 

30 The real part of the spatial domain bars 3 02 is then 

clamped in step 108 to square it up. In the frequency 
domain, this reintroduces some insignificant harmonics 400 
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(Fig. 4A) . The resulting square wave 402 in the spatial 
domain is seen in Fig. 4B, this is the calibration target. 

The amplitude of the calibration target is then 
scaled up by multiplication and rounding to a desired 
5 dynamic range, for example 0-255 for an 8-bit 

representation. The target is then cropped to the desired 
print size in the sample direction. In this example if 
the target is to occupy 10 inches at 3 00 dpi, then the 
entire generated target which is 3000 samples long is 
10 retained. 

Fig. 5 illustrates a lenticular printed image 500 
with the calibration target 402 printed 110 to the back of 
a lenticular sheet 502 . The target 402 is visible to the 
left eye 504 without optical banding artifacts, but is 

15 absent to the right eye 506. As a result accurate 

calibration can be seen to have been achieved, and the 
lenticular sheet can be seen to be produced with a pitch 
corresponding to 20.5 lines-per-inch. The printer can 
then be set up to print the image strips accordingly. 

20 Should optical banding artifacts have appeared in the 

calibration target this would indicate that the lenticular 
sheet was not produced to the calibrated pitch of 20.5 
lines-per-inch . 

Referring again to Fig. 1, the following MATLAB 6.0 

25 syntax pseudo-code is offered to show one way to implement 
process 100. 
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TABLE I 

Step 101 generates a square wave, 
factor = test_lines-per-inch/dpi ; 
This line determines the ratio of a test pitch, in lines- 
5 per-inch, of the sheet to the resolution, in dpi, of the 
printer . 

n = dpi*prec; 

This line determines the product of the printer resolution 
and precision of the calibration. 
10 t [0:n-l]; 

This line determines the samples. 

tmp = mod(2n *f actor*t , 2n) ; 
This line generates a square wave of required frequency. 
w 0 = 2n *0.5; 

15 This line computes the normalized frequency for breaking 
up the interval (0,2n). 

yl = (tmp>=w G ) ; 
This line assigns a value of '0' to normalized t between 
(0,w o ), x l' elsewhere so the square wave can be seen. 
20 A step 102 is a fast Fourier transform: 

Yl=fft(yl); 
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TABLE II 



A step 104 is the removal of aliasing frequencies by 
retaining only the fundamental and odd harmonics up to the 
Nyquist frequency: 

n_freq = n/2 ; 
This line determines the Nyquist frequency. 

n_period = round ( test_lines-per-inch*prec ) ; 

spike = n period; 
These lines determine the location of the fundamental. 

k = 0 ; %initialize harmonic counter 

Y=zeros (size ( Yl ) ) ; 

Y(1)=Y1(1) ; 
These lines create an array of zeros 

while ( (2*k+l ) *spike < n_freq) 

spike_Jharmonic = ( 2 *k+l ) *spike ; 

Y(spike_harmonic+l) =Yl ( spike__harmonic + l ) ; 

Y ( n - spike_harmoni c + 1 ) =Y1 (n- ( spike_harmonic ) +1) ; 

k = k + 1; 

end 

These lines copy into the array the DC component, 
fundamental and odd harmonics up to the Nyquist frequency 
and their conjugate locations to remove the aliased 
frequencies 



TABLE III 



A step 106 is the inverse Fourier transform: 
yy=real (if f t (Y) ) ; 

30 " ^ " " ~ ~' ~ " ~ ~" — 
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TABLE IV 



A step 108 is the removal of small oscillations from the 
inverse transform. This is done by setting a threshold 
and removing oscillation from the ideal 0 and 1 amplitudes 
(for a square wave) by clamping oscillations within the 
threshold to 0 and 1 respectively: 

a=find((cal > 1) | ( (cal<l) & (cal>0.9))) ; 
cal(a) = 1; 

These lines clamp the oscillations around 1 of plus or 
10 minus 0 . 1 to 1 

b=find((cal < 0) | ( (cal > 0) & (cal<0.1))) ; 
cal(b) = 0; 

These lines clamp the oscillations around 0 of plus or 
minus 0.1 to 1 . 
15 Once this one dimensional array is generated, it may be 
converted to strips by replicating the array to form bars 
or strips. 

In an alternative to step 100 the square wave is over 
sampled at a rate greater than the actual sampling rate (A 
20 dpi) of the system. E.g., the square wave is simply 

generated assuming a much greater B dpi (say M times the 
actual dpi) . 

Step 102 is dispensed with and a spatial domain 

25 filter is then designed for step 104, removing frequency 
components beyond A/ 2 dpi the Nyquist frequency of the 
actual system. Such a filter may be designed using 
conventional methods such as the Remez algorithm using 
software packages such as MATLAB. The Remez exchange 

30 algorithm is an iterative procedure consisting of two 

steps. One step is the determination of candidate filter 
coefficients h{n) from candidate "alternation 
frequencies, " which involves solving a set of linear 
equations. The other step is the determination of 

35 candidate alternation frequencies from the candidate 

filter coefficients. The algorithm converges very fast, 
and is widely used in practice to design optimal filters. 

The signal is then subsampled to the desired sampling 
rate by retaining every Mth sample and discarding the rest 

40 to arrive at the end of step 106. 
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This example performs the filtering in the spatial 
domain as opposed to the frequency domain. 

The post processing to remove oscillations step 108 
may be performed as explained for the first example. 
5 The description of the embodiments is provided to 

enable any person skilled in the art to make or use the 
present invention. While the invention has been described 
with respect to particular illustrated embodiments, 
various modifications to these embodiments will readily be 
10 apparent to those skilled in the art, and the generic 
principles defined herein may be applied to other 
embodiments without departing from the spirit or scope of 
the invention. The present embodiments should be 
considered in all respects as illustrative and not 
15 restrictive. Accordingly, the present invention is not 

intended to be limited to the embodiments described above, 
but is to be accorded the widest scope consistent with the 
principles and novel features disclosed herein. 

What is claimed is: 



